Analytic solutions of the geodesic equation in axially symmetric space— times 
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The complete sets of analytic solutions of the geodesic equation in Taub-NUT-(anti-)de Sitter, 
Kerr-(anti-)de Sitter and also in general Plebariski-Demiariski space-times without acceleration are 
presented. The solutions are given in terms of the Kleinian sigma functions. 
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INTRODUCTION AND MOTIVATION 

Supernovae and other data indicate the existence of 
Dark Energy which is appropriately described by a cos- 
mological constant. Galactic rotation curves, gravita- 
tional lensing and structure formation indicate the ex- 
istence of Dark Matter which in some approaches is ac- 
counted for by assuming a modified gravity theory. In 
any case, properties of the gravitational field, of the un- 
derlying geometry, can only be explored through the mo- 
tion of test objects like point particles and light rays. 
Therefore analytic solutions of the geodesic equation giv- 
ing the orbits of test objects in standard and generalized 
gravitational models is an important issue. 

Analytic solutions of the geodesic equation are also of 
use for black hole physics and the creation of gravita- 
tional waves: they provide an analytic determination of 
the shape of black holes with a link to the black hole 
parameters also in generalized scenarios. They can 
also be used for the calculation of gravitational wave tem- 
plates for Extreme Mass Ratio Inspirals (EMRI) for 
systematically finding homoclinic orbits |3(, also in gen- 
eralized scenarios, for analytic solutions and systematic 
studies of effective one-body equations of motion [4| , for 
testing numerical codes for binary systems, and finally 
may also be of use for practical applications like geodesy. 

Analytic solutions of the geodesic equation were only 
known for the motion in Schwarzschild [f|, Reissner- 
Nordstrom, Kerr, and Kerr-Newman space-times [f|. 
Based on new techniques developed in [3] we derived the 
complete set of analytical solutions of the geodesic equa- 
tion in Schwarzschild-(anti-)de Sitter space-time and ap- 
plied it also to the case of spherically symmetric space- 
times in higher dimensions [8| . Here we present complete 
sets of analytic solutions for a wide class of stationary 
axially symmetric gravitational fields. 

The metric of a general stationary axially symmetric 
black hole gravitational field in four dimensions has the 
form 



ds 2 



A, 



(dt - Ad<p)' 



dr' 



A n 2 
~ sin 2 0(odt - Bd^) 2 - 4-dtf 2 . 



where the functions p 2 , A r , A#, A, and B will be speci- 
fied later. We derive all analytic solutions of the geodesic 
equation in various space-times of this type, namely 
Taub-NUT-(anti-)de Sitter (TNdS), Kerr-(anti-)de Sit- 
ter (KdS) and even in Plebanski-Dcmianski (PD) black 
hole space-times without acceleration. The set of PD 
black-hole solutions exhausts all electrovac type D space- 
times. Since the Hamilton-Jacobi equation is separable 
only in the subset of space-times without acceleration [§] , 
we conclude that we can analytically solve the geodesic 
equation in all electrovac type-D space-times without ac- 
celeration. 

In all stationary axisymmetric space-times we have a 
conserved energy and angular momentum 
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-f(t- Acp) -a^- sin 2 •&( at- Btp) (2) 
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-A^(f-A<p) 
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B^- sin 2 d{ai - Bip) , (3) 
pi 



where the dot denotes the derivative with respect to s. 

Orbits in Taub-NUT-(anti-)de Sitter space-times 
The stationary axisymmetric TNdS space-time is given 
by 0,111 

p 2 :=r 2 +n 2 , A:=2ncos'd, B := p 2 , A# := 1 (4) 

A r := r 2 - 2Mr - - (r 2 + 3n 2 f - n 2 (l - 4Ati 2 ) , (5) 
3 

where A is the cosmological constant and M and n are 
the mass and NUT charge of the gravitating body. 

With dimcnsionless quantities f = r/rg (rg = 2Af), 
n = n/rs, L = L/rg, and the substitution £ = cosi? the 
Hamilton-Jacobi equation separates and reduces to 



df 
dr 

where 



R, 



dr 
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L - 2nE£ 



R = (f 2 - 

e e = ae 



h 2 ) 2 E 2 



A r {Sf 2 + L A + k) 



(6) 
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Here, k is the separation constant (Carter constant), a = 
—(k - Sn 2 + AE 2 n 2 + L 2 ), [3 = AEnL, ~/ = k-5n 2 , and 



A,. 



A r ro . Here and throughout the paper we use the 
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FIG. 1: TNdS: Parameter diagrams for n = 0.5, A = 2.9-KT 5 . 
The gray scales denote the number of positive real zeros of the 
polynomial R: black = 4, gray = 2, white = 0. Parameters 
in the dashed region are forbidden due to @. 



Mino-timc r defined through p 2 dr = r$ds 5 = 1 for 
time-like and 8 = for null geodesies. 

For given A and n, the orbits are characterized by E, 
L, and k. The range of these parameters is restricted 
by the condition that R and 0£ have to be positive or, 
equivalently, possess a certain number of real zeros. 

<d^ has two real zeros £1,2 € (— 1, 1) and, thus, can be 
positive only if 



ci 



Sh 2 + 4E 2 n 2 > 



c-2 = k — 5h 2 



L 2 



> 



(9) 



which constrains the allowed values for E and L. Then 
£ is restricted to the interval [^1,^2] which for n 5^ is 
not symmetric with respect to £ = 0. Thus, the motion 
of a test particle in TNdS space-time is bounded by two 
cones with in general different opening angles 12j. If 
k-8n 2 > then d g (0, tt) (see ng |2(b)| ), but for k-Sn 2 < 
the motion takes place in one half space (1? € (0, §) or 

g (f ,tt)) only (see Fig. [2(a)] ). If n = or L = 
the opening angles of the two cones are equal and the 
motion is symmetric with respect to the equatorial plane. 
It also can be shown along the lines of [12J that all test 
particles move on a cone which is tilted with respect to 
the symmetry axis and which touches the two bounding 
cones related to £1 and £2, see fig. [5] 

The function R is a polynomial of 6 th order in r which 
for positive A possesses at most 4 positive real zeros. The 
r coordinate may take negative values. 

Two examples of (E, L)-diagrams for n = 0.5 and dif- 
ferent k are shown in fig. [1] where the dashed region de- 
note forbidden values for E and L due to ©. Black 
indicates an escape orbit, a bound orbit (e.g. fig. Qtb)[ ) 
and a bound orbit crossing r = 0, gray one escape and 
a bound orbit crossing periodically r = 0, and white an 
escape orbit crossing r = (fig. [ ^a)[ ). 

With © we have a < and D := 1 - Aaj > and 
the solution for ■& is 



$(t) = arccosf — 



(VD sin (tc? - V^r) - /?)) , (10) 



(a) escape orbit crossing r = 0; (b) bound orbit; E 2 = 0.94038; 
E 2 = 0.99515, k = -0.5, k = 1.0, L = 2.0 

L = 2.5. 

FIG. 2: TNdS: Orbits for n = 0.5, A = 2.9 • 10~ 5 . 



where Tq = ^J—oltq + arcsin 2a ^°^P . A substitution f 

• + 77?, where ffj is a zero of i?, reduces the equation 
for f in (JSJ) to the hyperelliptic differential of first kind 

dr = f dx for a polynomial P5 of degree 5. This can 

V p 5(z) 

be solved by [7j 



(11) 



where cr^ are the derivatives of the Klcinian a function, 
/ describes the 0-divisor, i.e. a((f(z),z) t ) = 0, and 
Tn = Tn + f°° f , ■ The solution of the w equation is 



<p(r) 




arctan ■ 



1 - uB_ 



-A- arctan ■ 



1 - uB, 



where 



L±2Eh 



Bl - I, 



and B± 



B 2 
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f3±2a 





(12) 



' |L±2£nl"T 

For a detailed discussion see [15j. 
Orbits in Kerr-(anti-)de Sitter space-times The sta- 
tionary axially symmetric KdS metric written in Boyer- 



Lindquist coordinates t 

1 + ±a 2 A, emerges from ((T|) for [l3l . Il4| | 



\ 



p 2 := r 2 + a 2 cos 2 d, A# := 1 
A r := (1 - ±Ar 2 ) (r 2 + a 2 ) - 2Mr , 
^^asin 2 ??, B:=r 2 +a 2 . 



±Aa 2 cos 2 $ . 



where 

(13) 
(14) 
(15) 



The Hamilton-Jacobi equation again separates and 
gives with a = a/r$ and V = L/E, k = k/E 2 , 8 2 = 5/E 2 
and the substitution v = cos 2 •& (we exclude v = 1 corre- 
sponding to the coordinate singularities i? = 0, w) 



C ^ = EVR, ^ = E^Q h 
dr dr 



cl7 = ^ 



2D) 5 sin 2 d - V 



A,? sin 2 $ 



(16) 
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FIG. 3: Different types of geodesic motion. Here rs = 2, 
a = 0.4, k — 4, and A = 10~ 5 . The gray scales denote the 
number of positive real zeros of the polynomial R: black = 
5, gray = 3, white = 1. Parameters in the dashed region are 
forbidden due to Q v > 0. 



where 

R := X 2 (r 2 + a 2 - aV) 2 - A r (S 2 r 2 + k) (18) 
6„ := (1 - i/)A*(/s - 6 2 d 2 is) - x 2 (5(1 - v) - V) 2 (19) 

Solutions of HH) require Q u (v) > and Rif) > 0. 
Therefore, the number of real zeros of 0„ (in [0, 1)) and 
R depending on the parameters of the black hole and the 
particle gives all possible types of orbits. 

The number of zeros of 8„ in [0, 1) changes if (i) or 
1 is a zero or (ii) two zeros coincide. We have case (i) if 
L = aE±^ or L = 0, respectively. Considering case (ii) 
we can solve 0„ = (v — u) 2 (aiv + ao) with some constants 
a» for E 2 and L, where only u <G [0, 1) is of interest. In 
the same way we can analyze the changes in the number 
of zeros of R with the ansatz R= (r — u) 2 (X)i=o a i^)- 

A typical result of such an analysis for timclikc 
geodesies is shown in fig. [31 The dashed region denotes 
forbidden parameter values. The white area corresponds 
to an escape orbit, the gray area to an escape orbit but 
also a bound orbit, and the black area to an escape and 
two different bound orbits. 

The differential equation for $ is of elliptic type and 
first kind and is solved by 



$(r) = arccosi. 



di 



4p(2£V - f ;g 2 ,g 3 ) - d 2 



(20) 



(a) escape orbit, E 2 = 0.98, (b) bound orbit, E 2 = 0.94, 
L = X. L = 0.6. 

FIG. 4: KdS: Orbits for a = 0.4, A = 10~ 5 , and k = 3. 



and fn is a zero of R and 05 > 0. With f = E^/o^tq 
j°° the solution is given by 0] 



fM ^2 / / {Ey/a^T — fo ' 
[ ' + a 1 \ E^EEt-t 



The (^-motion depends on both r and 1? 

X~ 2 (tp - (fo) = I r - h , 



(22) 



(23) 



where we have to solve an elliptic and hyperclliptic inte- 
gral of third kind [l5| 



2ci 

R 



T * = ( S ^TT fc(^)(« - «o) + log - 



a(v - Vi) 



(6n + 5 i2 )a 3 A(x - aAV) + (6 i3 + S U )V 



(v - Vi), 
(~di) 



4x 



+ (a-V)(v-v ) 



(24) 



where v(t) = 2Et - f , v a = v(t ), p{v\) = \(a 2 Adi - 



d 2 ) 



(v 2 ), p(v 3 ) = z(di + d 2 ) = p(vi), and 



If = C 1 (w-w ) + C (f(w) + f(-ivo)) + J2 



05\XQ\ 

axo 



2.i 



Si{w) - Si(w ) - 




, (25) 



d 3 di, g 3 — \d\d 2 d 3 — jgS 2 a 4 d 2 — \d\, where w(t) = E^/a^r — fo, wo = w(tq), and Si(w) 



dy' 



and di = k — x ifl — ^ ) , 



where g 2 = \d\ 
f = 2Et + C 

V0 92J'~J3 

3d 2 = d 2 {nA -6 2 )-K + 2 X 2 d(d - V), Ad 3 = d 2 (6 2 {l - 
Aa 2 ) — x 2 — A/c) . Contrary to the TNdS-case, this motion 
is symmetric with respect to the equatorial plane. 

The differential equation for f is of hyperclliptic type 
and first kind. It can be cast in standard form x^- = 



log .((/M,^-2,r£ d z) _ Hcre d2 

CT ((/(«-),«,) t -2/^ dz) 

dr = {drudrrf with d n = El=i(k + 1 - *) a "+ 1+i ^ 
(see [IBl). The constants Cj are given by a partial frac- 
tion decomposition. The points pf = (xi , ± \/R Xi ) are 
the zeros of on the Riemann surface of y 2 = R x . 

Two orbits are shown in fig. [U for a detailed discussion 
see 0. 

Orbits in Plebariski-Demiariski space-times The PD 
space-times without acceleration are described by ([T]) 



Ey/a 5 R x by a substitution f = ±l/x + ffj where 
R x = f 



i=0 



a-. 



en 



(i"R) , (21) 



(6-i)! df( 6 " 1 ) 



4 



with [H 



p = r + (n — a cos -d) , 

= 1 + ±a 2 A cos 2 $ - | Aan cos $ 
A r =r 2 - 2Mr - n 2 + a 2 + Q\ + Q 2 n 

- |A (r 4 + (6n 2 + a 2 )r 2 + 3(a 2 - n 2 )n 2 ) 
A = a sin 2 # + 2n cos i? , B 



r 2 +a 2 +n 2 



(26) 
(27) 

(28) 
(29) 



where Q a and Q m are electric and magnetic charges of a 
gravitating source. Also in this case the Hamilton Jacobi 
equation separates (for simplicity, we choose neutral test 
particles; charged particles do not change the structure 
of the equations) and yields differential equations for f , 
£ = cos$, and <p 



df 

7w 

with 



d<^ 51 L — AE 



A,, sin 2 



(30) 



R = X 2 — A r (Sr 2 + k) (31) 
9 C = A*(l - i 2 ) (k - 6{n - 5£) 2 ) — (L— AE) 2 , (32) 



where X = (r 2 



a 2 



ah and A = Ar^ 



A standard substitution £ = ±i + £e> where £e is a 
zero of Qj, brings the equation for £ to a holomorphic 



hypcrelliptic differential of the first kind dr 
The solution of the equation for in (f3"0")) is 



ydy 



0(r) 



arccos 



o-2 / /(t - 



£e 



As in pip we obtain for the f motion the solution 



f(r) 



(33) 



(34) 



The function / describing the 0-divisor depends on the 
definition of the cr-function and, therefore, on the poly- 
nomials Ph{y) for the i? motion and on P^{x) for the f 
motion. 

<f depends both on i? and f giving a hypcrelliptic inte- 
gral of third kind which can be solved completely [IH . 

Therefore we succeeded in obtaining the complete an- 
alytic solution of the geodesic equation in all PD black 
hole space-times without acceleration. The orbits now 
depend on the three particle parameters E, L, and k 
as well as the six parameters characterizing this class 
of space-times. It has been shown that the PD black 
hole solutions exhaust all electrovac type D solutions. 
Since the condition of a space-time of being of electrovac 
type D without acceleration ensures separability of the 
Hamilton-Jacobi equation 0] we, thus, arrived at the 
conclusion that we now can explicitly give all analytic 



solutions of geodesic equations in all electrovac type D 
space-times without acceleration. 

In all cases discussed here the motion can be extended 
through r = to negative r [HI]. 

For bound orbits, the r and d motions are related to 
characteristic periods given by 



dr 



(35) 



which arc related to the zeros of the underlying polynomi- 
als. From these periods one can derive the perihelion shift 
and the Lense-Thirring effect, and also perform further 
orbital frequency analysis [17| . For escape orbits, the two 
deflection angles are given by Ai? = lim ($(s) — s)) 

s — >oo 

and Aip = lim (ip(s) — <p{— s)). These observable peri- 

s — >oo 

ods and deflection angles can be calculated analytically. 
For detailed calculations and discussions see [l5j |. 
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